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Abstract 

We first review the various definition of the total energy in the 
gravitational system. The naive definition has some defects, and we 
review how to modify the definition of the total energy. Then we 
explicitly demonstrate how to calculate the total energy of the system. 
Our example is the total energy of a black hole in the expanding 
closed de Sitter universe in (2+1) dimension. In general, we find that 
the contribution to the total energy comes only from the singularity. 
Then we can calculate the total energy by evaluating the contribution 
around the singularity. 



1 Introduction 

It has a long history how to define the total energy in the gravitational theory. 
Because of the invariance under the general coordinate transformation, it 
becomes ambiguous how to define the time variable. Then the Hamiltonian, 
which is conjugate to time variable, becomes ambiguous. Especially there 
are many ways||l|, 0, |^, |] how to define the energy of the gravitational field. 
The purpose of this paper give the prescription to calculate the total energy 
of the system with local singularity. It is known that (i) -Etotai = (finite) for 
the asymptotically fiat space-time and (ii) -Etotai = for the closed universe. 
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The contribution to the total energy comes only from the singularity, so 
we can calculate the total energy by evaluating the contribution around the 
singularity. 



2 Energy of Gravitational Field - Review - 

2.1 Action and the Einstein's Equation 

The Hilbert-Einstein action is given by |^ 

/ = ^ y" d'^Xyf^ R + j d^X^/^ >Cmatter , (l) 

where k = 8ttG. Taking the variation with respect to g^^, we have the 
Einstein's equation SIuE/^dfiu = 0, which gives 

R^"" - ^Rg^"" = kT^^ (2) 

with the boundary condition 

J d^x d^B^" = 0, (3) 

B'' = V^{g'"'iSTi,)-g-^\SK,)}. (4) 

Notice that the boundary term B'^ contains not only the variation of the 
metric 6g^i, but also the variation of the first derivatives of the metric Sg^^^x 
through the affine connection. 



2.2 Einstein prescription 

In the naive approach, the boundary term B'^ contains not only the variation 
of the metric but also the variation of the first derivatives of the metric, 
which makes difficult to find the classical solution compatible with boundary 
condition. There are various approaches to avoid this problem. 

First we review the primitive expression of energy conservations derived 
by Einstein in order to avoid the appearance of the variation of the deriva- 
tive of the metric in the boundary. Using the identity 

^R = ^G + d,V^^, (5) 



2 



G = ^^'^^{r^.rj, - r;;/eA}, (6) 

and the action is taken in the following form 

/* = ^ / d^xV^G + J d^xV^C^^tter. (8) 

After taking the variation of /* with respect to g^,^, the same Einstein's 
equation Eq.(^ is obtained with new boundary condition 

J d^xd^B*!" = 0, (9) 
B*^ = -V^g [5{V^ggnKx " 5{V^g''^)V^,^} . (10) 

In this expression, the boundary term B*^ contains only the variation of 
the metric as expected. Putting C*q = \J —gGj (2k), the energy- momentum 
tensor of the gravitational field is given by 

ciC* 

< = Q^g^ - ^'^h- (11) 

Then the conservation law of energy-momentum is expressed as 

a^« + T^') = 0, (12) 

where denotes the energy-momentum tensor density of matter. However 
the separation of the dynamical metric from the auxiliary metric is not clear 
in this method. 



2.3 ADM canonical formalism 

In order to clarify the physical meaning of the metrics, we often use the ADM 
formalism where we decompose the 4-dimensional space-time metric into 
(3+1) dimensional space and time metric. We consider the Cauchy surface 
Sf parametrized by a global function of time t. Let n"- be the unit normal 
vector to the hypersurface S^. The metrics are decomposed into the form 

ds^ = -N^dt^ + hab{dx'' + N''){dx^ + N''), (13) 
-N^ + N-Nt,NA , . 
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where is the lapse function, A^'^ is the shift operator and hab is the spatial 
metric. In this expression, the determinant of the metrics is g = Nh. 

First we proceed the Lagrangian formulation under the ADM (3+1) de- 
composition. Defining the extrinsic curvature of 

K,t = hy,nb (=ArrOj (15) 

and its mean value K = K^, the gravitational action is taken in the form 

Ig = ^{fd^xy/^R+fd^xVh2K), (16) 

= — f d^xNVhiKabK'^'' - +'^''^ R). (17) 

The notation (3) in the left upper indices denotes the corresponding three 
dimensional quantities. Taking the variation of Iq we have the Einstein's 
equation with the boundary condition, which contains only the variation of 
the first time derivatives of metrics. The surface term / dx^dtK/ k plays the 
role to make the action well defined in the Euclidean section [0. 

Next we proceed to the Hamiltonian formulation. The canonical momen- 
tum and the Hamiltonian density are expressed as 



TT 



dhab 



./h{K°-^ - Kh"^), (18) 



= NH + N^Pa + 2Da{h-^^^Nb7r^^). (19) 

In the above expression, the last term contribute only to an unimportant 
boundary term and will be dropped in the following. The Hamiltonian func- 
tion and momentum function are given by 

Pa = -2^rh{Dah'"^7i-'). (21) 

Because the time derivatives of A^ and N"- do not appear in the Lagrangian, 
they are auxiliary variables. Then we must put the Hamiltonian and the 
momentum constraints as 

H ^Q, (22) 
Pa ~ 0. (23) 
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The dynamical equations, which are derived by the canonical formalism, are 
given by 

Aab, (24) 
-Bab, (25) 
where 

Hg = j d^x{NH + N'^Pa). (26) 

While, if we directly take the variation of this ADM Hamiltonian 

5Ho = j rf=^x(A,fe57r"^ + B'^'dhab) - SC, (27) 

we obtain an extra surface contribution 6C. So the dynamical equation 
derived by the canonical formalism and the dynamical equation derived by 
the Lagrangian formalism give the different boundary condition. 

In order to avoid the additional surface term 6C, Regge-Teitelboim p 
modify the total Hamiltonian in such a way as the dynamical equation, de- 
rived by the canonical formalism, is equivalent to the equation derived by 
the variation of the total Hamiltonian. The explicit form of 6C is 

6C= j d^SiG'^'^NSh^.^k - NkSh,,) 

+ j d^Si{2Nk5T^^^ + {2N''ti^^ - N^Ti^^)5hjk}, (28) 

where 

Qijki ^ ^j^i/2(^j^ikyi ^ f^iij^jk _ 2h'ih^^). (29) 

If this surface term would vanish, Hq would be the correct Hamiltonian. 
However for the stationary asymptotically flat space-time [hij — 6ij ~ 1/t), 
the first surface term is different from zero, that is 

fd^SiG'^''%N6h,,,k - N^kSh,,) ~ 6 fd^Si{hu, - (30) 



ab 



71 



G 



5h 



ab 
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A new gravitational Hamiltonian is defined by 



H'^ = HG + E[h,,l (31) 
where additional contribution comes from the surface term by 

E[K^] = l^SiG'^'^'K.^u. (32) 



This is called ADM energy for the stationary asymptotically flat space-time. 
We get for the variation in this new gravitational Hamiltonian 

5H'^ = Jd'x{Aab67r''' + B^'6hab), (33) 

and the correct equations of motion Eq.(0)-Eq.(^) are recovered. 

The value of the energy becomes in the following. First, one has Hq = 
because of the constraint equations Eq.(pO|)-Eq. (pT]). For the surface con- 
tribution one has E[hij] = for the closed universe, because there is no 
boundary, and E[hij] = (finite) for the asymptotically flat space-time. For 
example, for the Schwarzschild metric, 

ds"" ~ -(1 - ^^)de + {5i, + 2Gm^)dx'dx\ 



one has well known result E[hij] = m & 10 1 



2.4 Comment on the surface term 

The surface terms in Eq.(|16D and Eq.(|32|) are important for the Lagragian 
formalism and the Hamiltonian formalism. Here we comment these surface 
terms are related to the surface term appeared in the identity Eq.(^-Eq.(|^). 
A. The surface term, which is the mean value of the scalar curvature in 
Eq. (p!6|) , plays a similar role to the time component of the surface term Eq.(|^) 
in Einstein's prescription in order to avoid the time derivative contribution 
in the surface term; 
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B. Using the three dimensional version of the identity Eq.(|^) -Eq.(0) 



^/h^'^R = Vh^'^ G + di^''>v\ (34) 

= ^'^G^'^'hi^^k, (35) 

where 

(3)Q^Jkl ^ Kl/2^j^ikyl ^ j^ilyk _ 2h'^h^l]^ (36) 

the additional energy Eq.(|32D contributed from the surface term is simply 
expressed as 

E[K,] = ^ I d'xd/^'^v'. (37) 

From the above observation the contribution to the total energy is inter- 
preted in the following way. The term y/—gR includes the singular terms 
coming from second derivative terms, and ^/—gG contains the less singular 
terms including only first derivative terms and the term dfjT>^ picks up the 
singular part as surface terms. 

Next we study the cases of neither closed space nor asymptotically flat 
space-time. 



3 Energy of System with Conical Singular- 
ity 

3.1 (2+1) dimensional conical singularity 



We consider the system of the point particle with mass niQ [|TT|. Static 
spherically symmetric metric is given by 

ds^ = -de + e'^(^)(dr2 + rHct)^). (38) 

Einstein equation becomes 

Acr(r) = mo5(r). 
7 



Then the solution is given as 

ds^ = -dt' + r-^^"*" {dr^ + r^d(P'), (39) 

where the ranges of variables are < r < oo and < < 27r. After making 
the change of variables 



l-4Gmo 



/{l-AGmo), 0=(l-4Gmo)0. (40) 



we get the flat metric 

ds^ = -df + + f2rf0^ (41) 
0<f<cx), O<0<27r-A, (42) 

where the deficit angle A = SyrGmo appears in the range of variable (j). This 
metric is not asymptotically fiat and the total energy formula Eq. (|3TD cannot 
apply. Here we define the total energy of this system as 

^totai = -7^ / d^xH, (43) 

where H is defined in Eq. (pO]) . Note that the minus sign, three dimensional 
integration and lack of the lapse function in -Etotai of Eq.(^). Using 
the explicit expression for Hin Eq. (PD|) , the correct value of total energy is 
obtained [pJ] 



^totai = 7^ / d^xVh^^^R =-^i dSiVa{x) = mo. (44) 

Next we explore the gravitational energy for the case with a singularity in 
the closed universe. 



3.2 Closed and expanding de Sitter universe with 
conical singularity in (2+1) dimension 

We consider the system of closed and expanding de Sitter universe with 
conical singularity. The classical solution is given by |J12| 

ds^ = -de + aitf (y3^ + (1 - mm)f''d(lA . (45) 
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Here the scale factor of the universe a{t) is determined by the equation 
a} jo} = —1/a^ + X , and the solution is given by a{t) = cosh{y/\t) / y/X. 
After making the change of variable = (1 — 6)0, the metrics in Eq.(^) 
shows the appearance of the deficit angle A = 27ib ~ SnGm, where the 
parameter 6 = 1 — yl — 8Gm is introduced. In order to evaluate the total 
energy of the system we make the change of variables 

f = r^-V(l -6), 

and we get the metric in the form 

d,' = -A^ + + r^d^^) . (46) 

The total gravitational energy is defined as the same form in Eq.(^). Almost 
all region in the integration vanish due to the constraint Eq.(^) and only 
contribution comes from the singular part around r = 0. Therefore the inte- 
gration is performed small region e around r = and the total gravitational 
energy is obtained as 

K Jo Jo or \ or ^ j 
= 2m, (47) 

where A is the cosmological constant. Note that the value of the total energy 
for the closed expanding de Sitter solution with conical singularity in Eq.(|4^) 
is twice for the simple conical singular solution in Eq. (^4D . 

The gravitational energy in asymptotically de Sitter spaces is examined 



by Abbott and Deser |13|. The three dimensional cosmological gravity is 



studied by Deser and Jackiw and found static many-body solutions . Mass 



in static asymptotically de Sitter spaces was studied by V. Balasubramanian 



et al. |]T5|. 



4 Summary and Discussion 

We first review the derivation of the total gravitational energy, i) Einstein 
prescription, ii) ADM canonical formalism using Regge-Teitelboim deriva- 



9 



tion. 

In section 3, we explicitly calculate the total gravitational energy for nei- 
ther closed space nor stationary asymptotically fiat space-time. Examples are 
the following (2+1) dimensional systems, i) system with conical singularity, 
ii) closed and expanding de Sitter universe with conical singularity. 

In general, we find that the contribution to the total energy comes only 
from the singularity. Then we can calculate the total energy by evaluating 
the contribution around the singularity. We will show the complete proof for 
our statement near future. 
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